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Abstract 

A binary supermassive black hole leaves an imprint on a galactic nucleus in the form of a "mass deficit," a 
decrease in the mass of the nucleus due to ejection of stars by the binary. The magnitude of the mass deficit 
is in principle related to the galaxy's merger history, but the relation has never been quantified. Here, high- 
accuracy A^-body simulations are used to calibrate this relation. Mass deficits are shown to be M^igf w 0.5Mi2, 
with Mi2 the total mass of the binary; the coefficient in this relation depends only weakly on M2/M1 or on 
the galaxy's pre-existing density profile. Hence, after 9V; mergers, Mi/gj w 0.59V^M, with M, the final (current) 
black hole mass. When compared with observed mass deficits, this result implies I ^ 9\[ <; 3, in accord with 
hierarchical galaxy formation models. Implications for binary stalling radii, the origin of hyper-velocity stars, 
and the distribution of dark matter at the centers of galaxies are discussed. 
Subject headings: 



1. INTRODUCTION 

Galaxy mergers bring supermassive black holes (SBHs) 
together (iBe^elman. Blandford. & Rees ' 1980), and binary 
SBHs are increasingly invoked to explain the properties 
of normal and active galaxie s (Komossa 2003), including 
AGN variability ('Valtaoja et al. 2000; Xie 2003), the bending 
and p recession of radio jets (Roos et al. 1993; Romero et al. 
20001). X- and Z -shaped radio lobes ( Merritt & E kers 2002; 



Gopal-Krishna et al. 2003) and the presence of cores in brigh t 



elliptical galaxies dMilosavlievic et al.ll2002l; iGrahamI 12004 . 
Binary SBHs may also be responsible for t he high-velocity 
stars observed in the halo of the Milky Way (lYu & Tremaind 
l2003t iHaardt et al.l l2006'). and for other populations that ap- 
pear to have been ejected fr om galaxies, including int ra- 
cluster planetary nebulae (Hol ley-Bockelmann et aljr2005h . 

The dynamical interaction of a massive binary with stars in 
a galactic nucleus is often discussed with refer ence to rate 
equations d erived from scattering experiments ( HillsllT983t 
iMikkola &^ V altonen 1992; Quinlan 1996.). A massive binary 
hardens at a rate 



dt \ a J 



(1) 



where a is the binary semi-major axis, p and o are the stellar 
density and velocity dispersion, and // is a dimensionless rate 
coefficient that depends on the binary mass ratio, eccentricity, 
and hardness. The mass in stars ejected by the binary satisfies 
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where M12 = Mi +M2, the binary mass, and 7 is a second 
dimensionless coefficient. Equations ([T]i and (|2|i have been 
used to estimate evolution rates of binary SBHs in galactic 
nuclei (e.g. Volonteri et al. 2003a, b; Sesana et al. 2004) 
and to compute the mass in stars ejected from the galaxy by 
the binary (e.g. Holey-Bockelmann et al. 2005; Haardt et 
al. 2006). These equations are also the basis for many hybrid 
schemes wh ich imbed a binary into a model of the host galaxy 
jZier & Bie rmann 2001; Yu 2002; Merritt & Wang 2005). 

Equations ([TJ and ^ are not particularly useful however 
when computing the binary-induced change in the distribu- 
tion of stars in a nucleus, which is the topic of the current 



paper. A^-body experiments show that the two SBHs substan- 
tially change the stellar density in a short time, of order the 
galaxy crossing time or less, before they form a tightly-bound 
pair (e.g. Figure 4 of Milosavljevic & Merritt 2001). In other 
words, by the time that equations like ([T]i and (|2]) start to be- 
come valid, a considerable change has already taken place in 
the stellar distribution. Furthermore, in many galaxies, evo- 
lution of the binary would be e xpected to stall at about the 
same time that it becomes hard (IValtonenlll996h . Even if the 
mass ejected by the binary before this time could be accu- 
rately computed from an equation like (|2]), it would still be 
difficult to relate Mej to observable changes in the nuclear 
density profile. 

A^-body techniques would seem to be the solution to these 
problems; but A^-body integrations are plagued by spurious 
relaxation and other discreteness effects, which cause an em- 
bedded binary to continue evo lving even in circumstances 
where a real binary would stall jBerczik. Merritt & SpurzemI 
l2005h . The counterpart of these discreteness effects in real 
galaxies is two-body relaxation, but - especially in the lumi- 
nous galaxies that show evidence of the "scouring" effects of a 
massive binary - central relaxation times are much too long t o 
significantly affect the supply of stars to the binary (lYul2002h . 
A^-body simulations that are dominated by discreteness effects 
can not be scaled to the essentially collisionless regime of real 
galaxies. 

However, there has been much progress in the art of A^-body 
simulation in recent years, such that high-accuracy, direct- 
summation integrations are now feasible with particle num- 
bers as large as 10^ (Dorband et al.' 2003; Gualandris ( 
20041; pukushige et al. 2005; Berczik, Merritt & Spurzer 
2005h . Given such high values of A^, it becomes possible to 
separate the rapid, early phases of binary evolution - which, 
when appropriately normalized, should be independent of A^ 
- from the late, slow stages that are driven by A^-dependent 
process like collisional loss-cone repopulation. 

This is the approach adopted in the present paper. A^-body 
simulations are used to follow the evolution of a galaxy con- 
taining a central massive object and a second, inspiralling 
point mass. The integrations are continued until the two mas- 
sive particles form a tight binary. The time at which the evo- 
lution of the binary switches from rapid - i.e. A^-independent 
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- to slow - i.e. A^-dependent - is identified, and the properties 
of the galaxy are recorded at that time. The degree to which a 
clear separation of the two regimes is possible, for a given A^, 
depends on the binary mass ratio, becoming more difficult as 
the mass ratio becomes more extreme. Results are presented 
here for mass ratios in the range 0.025 < MijMx < 0.5. 

We find (§4) that the mass deficit - the difference in inte- 
grated mass between initial and final nuclear density profiles 

- at the end of the rapid evolutionary phase is proportional to 
the total mass of the binary, M^i^f « 0.5Mi2, with only a weak 
dependence on M2/M1 or on the initial density profile. This 
result can be motivated by simple arguments (§2): the smaller 
M2, the tighter the binary which it forms before stalling. 

A mass deficit of ~ 0.5M, is at least a factor two 
smaller than the typical mass de ficits observed in bright 
elliptical galax ies (Milosavlievic et al. 2002; Graham 2004; 
iFerrarese et al.l 12006) but we show via an additional set of 
A^-body simulations (§5) that the effect of binary SBHs on 
a nucleus is cumulative, scaling roughly in proportion both 
to the number of mergers and to the final mass of the SBH. 
Hence, observed values of M^/gf/M, can be used to constrain 
the number of mergers that have taken place since the era 
at which the SBHs first formed. We find (§6) that l<,9{, ^3, 
consistent with expectations from hierarchical structure for- 
mation theory. Finally, in §7 and §8, the implications for 
binary stalling radii, ejection of high-velocity stars, and the 
distribution of dark matter are discussed. 

2. STAGES OF BINARY EVOLUTION 

Here we review the stages of binary SBH evolution and dis- 
cuss the connection between evolution of a binary in A^-body 
simulations and in real galaxies. 

Let Ml and M2 to be the masses of the two components of 
the binary, and write Mn = Mi +M2, q = M2/M1 < 1, and 
fj = M1M2/M12. In what follows we assume that the larger 
SBH is located initially at the center of the galaxy and that the 
smaller SBH spirals in. 

The evolution of the massive binary is customarily divided 
into three phases. 

1. At early times, the orbit of the smaller SBH decays due 
to dynamical friction from the stars. This phase ends when 
the separation R12 between the two SBHs is ~ r/,, the gravita- 
tional influence radius of the larger hole. We define r/, in the 
usual way as the radius of a sphere around Mi that encloses a 
stellar mass of 2Mi : 



M*(r,,)=2Mi. 



(3) 



2. When R12 falls below ~ r/,, the two SBHs form a 
bound pair The separation between the two SBHs drops 
rapidly in this phase, due both to dynamical friction act- 
ing on M2, and later to eje ction of stars by the binary 
(iMilosavljevic & Merrittll20011) . The motion of the smaller 
SBH around the larger is approximately Keplerian in this 
phase; we denote the semi-major axis by a and the eccen- 
tricity by e. The relative velocity of the two SBHs for e = 



GM12 



and the binary's energy is 



GM1M2 
2a 



(4) 



(5) 



3. The rapid phase of binary evolution comes to an end 
when the binary's binding energy reaches ^ M^CT^, i.e. when 



a w fl/„ where a/, is the semi-major axis of a "hard" binary, 
sometimes defined as 



ah 



4o2' 



(6) 



If one writes r/, = GMi /o^ - equivalent to the definition (O 
in a nucleus with p 



ah 



- then ah can also be written 

q Hi 

(1+^)2 4- 
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The exact definition of a "hard" binary varies from author to 
author, and this vagueness reflects the difficulty of relating the 
evolution of an isolated binary to one embedded in a galaxy. 
An isolated binary in a fixed background begins to harden at 
an appr oximately con stant rate, {d/dt){l /a) w const., when 
a <, ait (IOuinlanlll996.) . But in a real galaxy, a hard binary 
would efficiently eject all stars on intersecting orbits, and its 
hardening rate would suddenly drop at a « a/,. This effect 
has bee n observed in A^-body experiments with sufficiently 
large Af (iMakino & Funatol2004l:rBerczik. Merritt & SpurzemI 
l2005h . 

Supposing that the binary "stalls" at a « a/„ it will have 
given up an energy 
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to the stars in the nucleus. In other words, the energy trans- 
ferred from the binary to the stars is roughly proportional to 
the combined mass of the two SBHs. This result suggests that 
the binary will displace a mass in stars of order its own mass, 
independent of the mass of the infalling hole. This prediction 
is verified in the A^-body simulations presented below. 

Since the "hard" binary separation a/, is ill-defined for bi- 
naries embedded in galaxies, we define here a more useful 
quantity. The stalling radius astaii is defined as the separation 
at which evolution of the binary would halt, d{\/a)/dt = Q,\n 
the absence of any mechanism to re-supply the stellar orbits. 
For instance, if the galaxy potential were completely smooth, 
the binary's decay would halt once all the stars on orbits in- 
tersecting the binary had been ejected or otherwise (e.g. due 
to shrinkage of the binary) stopped interacting with it. The 
motivation for this definition is the very long time scales, in 
luminous elliptical galaxies, for orbital re-population by two- 
body encounters. Even in galaxies with relatively short cen- 
tral relaxation times, the binary's hardening rate would drop 
drastically at a = « and the likelihood of finding the binary 
at a separation near astaii would be high. 

This definition implies a dependence of astaii on the mass 
ratio of the binary, but (unlike Equation |6ll it also implies a 
dependence on the initial density profile and shape (spherical, 
axisymmetric, triaxial) of the galaxy, since these factors may 
influence the mass in stars that can interact with the binary. 
The definition is operational: it can only be applied by "doing 
the experiment," i.e. imbedding the binary in a galaxy model, 
turning off gravitational perturbations (aside from those due 
to the binary itself), and observing when the decay halts. 

In a numerical simulation with finite A^, gravitational en- 
counters will continue supplying stars to the binary at rates 
much higher than those in real galaxies. One can still esti- 
mate astaii if the particle number is large enough that a clear 
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TABLE 1 

Parameters of the A'-body integrations 



y 


M2/M1 






"stall 




Qstall/'-; 




Mdef/Mi2 


0.5 


0.5 


0.264 


0.37 


1.49 X 10" 


-2 


4.08 X 10" 


-2 


0.60 


0.5 


0.25 


0.264 


0.32 


1.18 X 10" 


-2 


3.69 X 10- 


-2 


0.46 


0.5 


0.1 


0.264 


0.30 


5.38 X 10" 


-3 


1.82 X 10" 


-2 


0.33 


0.5 


0.05 


0.264 


0.29 


3.60 X 10- 


-3 


1.26 X 10- 


-2 


0.27 


0.5 


0.025 


0.264 


0.28 


1.91 X 10" 


-3 


6.82 X 10- 


-3 


0.24 


1.0 


0.5 


0.165 


0.24 


7.75 X 10" 


-3 


3.23 X 10- 


-2 


0.62 


1.0 


0.25 


0.165 


0.22 


4.83 X 10" 


-3 


2.25 X 10- 


-2 


0.54 


1.0 


0.1 


0.165 


0.19 


2.37 X 10" 


-3 


1.23 X 10- 


-2 


0.45 


1.0 


0.05 


0.165 


0.19 


1.36 X 10" 


-3 


7.27 X 10- 


-3 


0.39 


1.0 


0.025 


0.165 


0.18 


6.21 X 10" 


-4 


3.36 X 10- 


-3 


0.34 


1.5 


0.5 


0.0795 


0.13 


3.66 X 10" 


-3 


2.77 X 10- 


-2 


0.63 


1.5 


0.25 


0.0795 


0.11 


2.18 X 10" 


-3 


1.95 X 10- 


-2 


0.56 


1.5 


0.1 


0.0795 


0.10 


1.26 X 10- 


-3 


1.29 X 10- 


-2 


0.41 


1.5 


0.05 


0.0795 


0.10 


7.78 X 10" 


-4 


8.19 X 10- 


-3 


0.38 


1.5 


0.025 


0.0795 


0.09 


4.60 X 10" 


-4 


4.84 X 10- 


-3 


0.39 



change takes place in the hardening rate at some value of a. 
Better still, by repeating the experiment with different values 
of A^, one can hope to show that the rapid evolutionary phase 
comes to an end at a well-defined time and that subsequent 
evolution occurs at a rate that is a decreasing function of A^. 

One advantage of defining fl.s,a// in this operational way is 
that it eliminates the need for many of the qualitative distinc- 
tions that have been made in the past between the different 
regimes of binary evolution, e.g. "soft binaries" vs. "hard bi- 
naries," hardening via "dynamical friction" vs. hardening via 
"loss-cone draining," etc. ( Yu 2002). It also allows for mech- 
anisms that are not reproduced in the scattering experiments 
on which E quations (1) and (2) are based , e.g. the "secondary 
sUngshot" ( iMilosavlievic & Merit3l2003b . or the effect of the 
changing nuclear density on the rate of supply of stars to the 
binary. 

It would seem natural to compute Ostaii using so-called 
"collisionless" A^-body codes that approxim ate the gravita- 
tional potential via smoo th basis functions ('Clutton-Broc^ 
l97llvan AlbadaX van G orkom 1977; Hernquist & OstrikeiJ 



1992). Such algorithms have in fact been applied to the binary 
SBH problem ( , Quinlan & Hernquist 1997; Hemsendorf et al.i 
I200I IChatteriee et al.ll2003ir but with results that are not al- 
ways consistent with those of direct-summation codes or with 
the predictions of loss-cone theory (compare Chatterjee, Hen- 
quist & Loeb 2003 with Berczik, Memtt & Spurzem 2005). 
These inconsistencies may be due to difficulties associated 
with incorporating a binary into a collisionless code, ambigui- 
ties about the best choice of origin for the potential expansion, 
etc. Until these issues can be resolved, direct-summation A^- 
body codes seem a safer choice. 

3. MODELS AND METHODS 

The galaxy models and A^-body techniques used here are 
similar to those described in Merritt & Szell (2005) (hereafter 
Paper I). In brief, Monte-Carlo realizations of steady-state, 
spherical galaxy models were constructed using Dehnen's 
(1993) density law. The models contained an additional, cen- 
tral point mass representing the more massive of the two 
SBHs. The mass of this particle. Mi, was always set to 0.01 in 
units where the total mass in stars Mf,ai was one. (The Dehnen 
scale length ro and the gravitational constant G are also set to 
unity in what follows.) Stellar positions and velocities were 



generated from the unique isotropic phase-space distribution 
function that reproduces Dehnen's p(r) in the combined grav- 
itational potential of the stars and the central point mass. The 
second SBH was introduced into this model at f = 0, and 
given a velocity roughly 1/2 of the circular orbital velocity. 
In the models with y = 0.5 and y — 1.0, the initial separa- 
tion of the two massive particles was 1 .6, while in the mod- 
els with Y= 1.5 the initial separation was 0.5. Each realiza- 
tion was the n integrated forward using the A^ -body integrator 
described in Merritt, Mikkola & Szell (2005), an adaptation 
of NBODYl 'Aai-seth (1999) to the GRAPE-6 special-purpose 
computer Close encounters between the two SBH particles, 
and between the SBHs and stars, were regularized using the 
algor ithm of Mikkola and Aarseth (Mik kola & Aa rseth 199(J 
[19931) . 

Three parameters suffice to define the initial models: the 
binary mass ratio q, the central density slope of the Dehnen 
model Y = — lilogp/dlogr, and the number of particles A^. 
Five values of the binary mass ratio were considered: q = 
(0.5,0.25,0.1,0.05,0.025). For each choice of q, three val- 
ues of Y were used: 0.5, 1.0, and 1.5. Finally, each of these 
15 initial models was integrated with two different values of 
A^: 1.2 X 10^ and 2.0 x 10^. Mass deficits were computed as 
described in Paper I. 

The parameters of the A^-body integrations are summarized 
in Table 1 . Columns three and four give two different esti- 
mates of the binary's influence radius. The first, r/,, is the 
radius containing a mass in stars equal to twice Mi in the 
initial model. The second is the radius containing a mass 
in stars equal to twice Mi +M2 at a time t = tstaih the esti- 
mated stalling time. The second definition, r^'^, is the relevant 
one when comparing the results of the A^-body integrations 
to a real galaxy which has already experienced the scouring 
effects of a binary SBH. The other columns in Table 1 give 
estimates of the stalling radius flj,„// = a{t stall) and the mass 
deficit M^gy^(rs,a/;), derived as described below. 



4. RESULTS 

Figure [T] illustrates the evolution of the relative orbit in the 
five integrations with Y = 1 -5 and A^ = 2QQK. The three phases 
discussed above are evident. (1) The separation R12 between 
the two massive particles gradually decreases as dynamical 
friction extracts angular momentum from the orbit of M2. (2) 
When R12 ~ rij, the rapid phase of binary evolution begins. 
(3) When the separation drops to ^ a/,, evolution again slows. 
Figure [T] shows clearly that the separation of the two massive 
particles at the start of the final phase is smaller for smaller 
M2. 

Figure [T] suggests that the stalling radius is of order a/,, but 
better estimates of can be made by comparing integra- 
tions of the same model carried out with different A^. Fig- 
ure |2] makes the comparison for the model with Y = 1 and 
^ = 0.1. After 15 orbits, the mean separation between the 
two massive particles has dropped from ^ 1 to ^ 0.2 « r/, and 
the second, rapid phase of binary hardening begins. Prior to 
this time, the two integrations with different A^ find essentially 
identical evolution of Rn- As discussed above, the evolution 
during this phase is due to some combination of dynamical 
friction acting on M2 and ejection of stars by the increasingly- 
hard binary, both of which are A^-independent processes. The 
contribution of dynamical friction to the evolution can be es- 
timated using Chandrasekhar's (1943) expression for the dy- 
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Fig. 1. — Evolution of the binary separation in five W-body integra- 
tions with Y= 1.5 and A' = 200K; binary mass ratios are, from left to right, 
0.5,0.25,0.1,0.05,0.025. Vertical fines show the times identified as f = fj,„;;. 
The upper horizontal line indicates r/,, the influence radius of the more mas- 
sive hole in the initial model. Lower horizontal lines show c/;, as defined in 
Equation Q. The rapid phase of decay continues until a o;,, with the result 
that the binary's binding energy at the end of this phase is neai'ly independent 
of Ml- 

namical friction force: 

(Av||) = -4ni4nG^M2p) 1^ dvf V/(v/)//i(v,vy()?a) 

^i = ln^?T^' (9b) 
10 ifv<v/. ^ ' 

This is the standard approximation in which InA is "taken 
out of the integral"; v is the velocity of the massive object, 
V/ is the velocity of a field star, and //(v/) is the field-star 
velocity distribution, normalized to unit total number and as- 
sumed fixed in time. The result of integrating equation (|9]l, 
with InA = 5.7, is shown in the upper panel of Figure |2] as 
the red (thin) line. Chandrasekhar's formula accurately repro- 
duces the evolution of the relative orbit until a time f « 80, af- 
ter which it predicts that the separation should drop to zero at a 
finite time. However as shown in the lower panel of Figure |2] 
by f = 80 - roughly at the start of the second evolutionary 
phase - the infall of M2 has begun to displace stars and lower 
the central density, and the dynamical friction force must drop 
below the value predicted by equation (|9]l with fixed p and //. 

The inset in Figure |2] shows the inverse semi-major axis of 
the binary, I /a, versus time in the two integrations of this 
model with different A^. The two curves are coincident un- 
til a Uf, but diverge at later times. T his is the expected 
result (iBerczik. Merritt & Spurzemll2005h : once the hard bi- 
nary has ejected most of the stars on intersecting orbits, con- 
tinued hardening requires orbital repopulation which occurs 
on a time scale that is roughly proportional to A^. The dif- 
ference between the two integrations is also apparent in the 
lower panel of Figure|2] which shows mass deficits in the two 
integrations. Based on this comparison, the stalling radius - 
the value of a at which the evolution ceases to be independent 
of particle number-is Qstaii =a{t « 110) Ri2 x 10^^, and the 
mass deficit when a = Ostaii is M^^f w 0.5Mi. 
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Fig. 2. — Upper panel: Thick (black) line shows evolution of the binary 
separation in the W-body integration with Y= 1.0, M2/M1 = 0.1, W = 200K. 
Thin (red) line is the evolution predicted by the dynamical friction equation 
(9) assuming a fixed galaxy. Horizontal lines indicate r;, and a;,, the latter as 
defined in Equation Q. The inset shows the evolution of the inverse semi- 
major axis of the binary for this integration, and for a second integration with 
N =\ 20K. Lower panel: The mass deficit, as defined in the text, for the same 
two A'-body integrations. Lines show least-squares fits to the time intervals 
80 < f < 110 and 110 < f < 200. 

In order to make still more accurate estimates of fj,„// and 
(istuih the hardening rate. 



s{t) 



1 

dt \a 



(10) 



was estimated from the A^-body data by fitting smoothing 
splines to the measured values of a^' vs. t and differenti- 
ating. In a collisionless galaxy, s{t) would reach a peak value 
during the rapid phase of binary evolution, then drop rapidly 
to zero as the stars on intersecting orbits are ejected by the 
binary; tstaii would be the time at which i w 0. In the A^-body 
integrations, s will never fall completely to zero since gravi- 
tational encounters continue to scatter stars into the binary's 
loss cone. 

Figure |3] shows s{t) as extracted from the A^-body integra- 
tions with Y = 0.5. In each case, s reaches a peak value during 
the rapid phase of evolution and then declines. In the mod- 
els with largest M2, loss-cone repopulation is least efficient, 
and the hardening rates in the two integrations with different 
A^ "track" each other well past the peak. As M2 is decreased, 
loss cone repopulation becomes more efficient and the two 
s{t) curves deviate from each other at progressively earUer 
times with respect to the peak. 
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Fig. 3. — Evolution of the binary liardening rate .v = (d / dt)( \ / a) in N- 
body integrations with y= 0.5; binary mass ratios are, from left to right, 
0.5,0.25,0.1,0.05,0.025. Solid curves are for N = 200K and dashed curves 
are for N = 120K. Vertical solid lines are the estimates of ?,,„;;. Nearly- 
vertical dashed lines show the approximate, asymptotic behavior of ,s(/) in a 
purely collisionless galaxy, assuming l[ah) = (40,50,60). 

Even in the absence of loss cone repopulation, a binary in a 
collisionless galaxy would continue to evolve at late times as 
stars with progressively longer orbital periods reach pericen- 
ter and interact with it. The asymptotic behavior of s{t) in a 
spherical non-evolving galaxy in the absence of encounters is 
roughly 



s[t) 



1 

dt \a 



16%^K I f„,{E)dE (11) 



where fm{E) is the phase-space mass density of stars and K 
is a constant that defines the mean, dimensionless change of 
energy of the binary in one interaction with a star ( Yu 2002'). 
Equation (fTTT i equates the rate of change of the binary's en- 
ergy with the rate at which stars, moving along their unper- 
turbed orbits, enter the binary's influence sphere and extract 
energy from it. Depopulation of the orbits as stars are ejected 
is represented by progressively restricting the range of the en- 
ergy integral; the lowest allowed energy at any time t is that 
correspondi ng to an orbit with radial period P{E) — t. Fol- 
lowing 1^ (120021) . we set = 1 .6, the approximate value 
for a "hard" binary, and assume that the integration starts at 
a = a/,. Figure|3]compares the solution to Equation (fTTl i with 
the A^-body hardening rates for (y = 0.5,^ = 0.5). Clearly, the 
"draining" of long-period orbits contributes only minimally 
to the late evolution seen in the A^-body integrations, i.e. the 
evolution in phase three is driven almost entirely by loss-cone 
repopulation. 

Plots like those in Figure |3] were used to estimate t stall for 
each (y,M2/Mi). Unavoidably, these estimates are somewhat 
subjective. For the larger values of M2, tstaii was taken to be 
the time when the 5(f) curves for the two different A^-values 
begin to separate. For the integrations with smaller M2, col- 
lisional effects are more important, and the two s{t) curves 
separate even before the peak value is reached. In these inte- 
grations, t stall was taken to be the time at which in the = 
200K integration reached its first minimum after the peak. 
These estimated tstaii values are likely to be systematically 
larger than the true values in a collisionless galaxy. However 
the uncertainties in the estimated values of a stall = '^{tstaii) and 
Mdefih-taii) are probably modest, at least for the larger M2 val- 
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Fig. 4. — Stalling radius (a) and mass deficit at f = (b) as functions of 
the binary mass ratio. Curve in (a) is Equation (12); curve in (b) is Equation 
(13). 

ues. For instance. Figure |2] shows that the mass deficit in the 
= 200K integration of the (y = 1.0,^ = 0.1) model varies 
only between 0.50 Mjef/Mi < 0.65 for 110 < f < 140. In- 
tegrations with larger A^ will ultimately be required in order to 
improve on these estimates. 

The results are presented in Table 1 and Figure|4] The upper 
panel of Figure |4] shows Osiaii as a fraction of r^^, the binary's 
influence radius (the radius containing a mass in stars equal to 
Ml +M2) at tstaii- The measured points are compared with 



Ostall 



0.2- 



q 



(12) 



this functional form was motivated by Equation Q. Given 
the likely uncertainties in the estimated Ostaii values, Figure|4] 
suggests that there is no significant dependence of astaii/r]^ 
on the initial nuclear density profile. This result will be used 
below to estimate stalling radii in observed galaxies. 

The lower panel of Figure|4]shows mass deficits at f = tstaii 
as a fraction of M12. The line in this figure is 



M, 



def 



Mn 



: 0.70^' 



,0.2 



(13) 



This relation is a good fit to the y = 1 .0 and y = 1 .5 models, 
though it overestimates M^/g/ for y — 0.5. In any case, the 
dependence of Mdef on binary mass ratio is weak, consistent 
with the prediction made above. 

In galaxies with initially steep nuclear density profiles, 
p ^ r^^, 1.0 ^ y^ 1.5, Figure m suggests that mass deficits 
generated by "stalled" binaries should lie in the narrow range 



0.4;$ 



M, 



def 



Ml +M2 



j$0.6, 0.05^^;$ 0.5. 



(14) 
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Fig. 5. — Density profiles at / = fiM/;- Thickness of curves decreases from q 

To a good approximation, M^e/ « 0.5M,. 

Density profiles att~ tstaii are shown in Figure |5] It is in- 
teresting that none of the profiles exhibits the very flat, nearly 
constant density cores seen in some b right elhptical gala xies 
( lKormendvll985l:lLauer et~ani2002l:lFerrarese et alj|2006h . 

5. MULTI-STAGE MERGERS 

The weak dependence of Mdef/Mn on initial density pro- 
file and on q found above has an obvious implication: in re- 
peated mergers, mass deficits should increase cumulatively, 
even if expressed as a multiple of M,, the combined mass 
of the two SBHs at the end of each merger. If the stellar 
mass displaced in a single merger is ^ 0.5Mi2, then (assum- 
ing that the two SBHs always coalesce before the next SBH 
falls in) the mass deficit following 9\[ mergers with M2 ^ Mi 
is ~ 0.59\[M, with M, the accumulated mass of the SBH. 

Equation (fTsT i could be used to derive a more precise pre- 
diction for the dependence of the mass deficit on . But di- 
rect simulation is a better approach. To this end, multi-stage 
A^-body integrations were carried out. The starting point for 
each set of integrations was one of the y = 1.5, Mi = 0.01 
A^-body models described above, extracted at f = t^iaii- The 
two massive particles were replaced by a single particle with 
mass equal to their combined mass, and with position and 
velocity equal to the center-of-mass values for the binary. 
A second massive particle, with mass equal to the original 
value of M2, was then added and the model integrated for- 
ward until the new stalling radius was reached. The process 
was then repeated. In each set of experiments, the mass M2 
of the second SBH particle was kept fixed; in other words, 
M,, the accumulated central mass, increased linearly with 
the number of mergers, while the binary mass ratio de- 
creased with Three different values of M2 were tried: 
M2 = (0.005,0.0025,0.001). 

Figure |6] and Table 2 give the results. The cumulative ef- 
fect of multiple mergers is clear from the figure: Mjef/M, 
increases roughly linearly with for the first few mergers, 
reaching values of ~ 1 for 94^—2 and ^1.5 for ?v; = 3. For 
larger ?v; and for smaller M2, the increase of Mjef / M, with 
begins to drop below a linear relation. That this should be so is 
clear from the results of the single-stage mergers (Figure |4]i. 
Nevertheless, Figure |6] verifies that (at least for ^ ^ 0.1) hi- 
erarchical mergers produce mass deficits that scale roughly 



r r 

-0.5\a q = 0.025. Dotted lines are tlie initial density profiles. 



TABLE 2 
Multi-StageA'-body integrations 



Ml 


3v; 


M. 








0.0050 


1 


0.0150 


0.0095 


0.63 


0.63 




2 


0.0200 


0.024 


1.20 


0.60 




3 


0.0250 


0.041 


1.64 


0.55 




4 


0.0300 


0.057 


1.90 


0.48 


0.0025 


1 


0.0125 


0.0070 


0.56 


0.56 




2 


0.0150 


0.017 


1.13 


0.57 




3 


0.0175 


0.027 


1.54 


0.51 




4 


0.0200 


0.035 


1.75 


0.44 


0.0010 


1 


0.0110 


0.004 


0.41 


0.41 




2 


0.0120 


0.010 


0.79 


0.40 




3 


0.0130 


0.014 


1.10 


0.37 




4 


0.0140 


0.018 


1.29 


0.32 



linearly with both 9^ and M, for small 9^ . 

Of course, these results are only meaningful under the as- 
sumption that the binary manages to coalesce between succes- 
sive merger events. Infall of a SBH into a nucleus containing 
an uncoalesced binary would almost certainly result in larger 
values of M^^f jM, than those given in Table 2, for reasons 
discussed below. 

6. COMPARISON WITH OBSERVED MASS DEFICITS 

Mass deficits, computed from observed luminosity pro- 
files, have been pubhshed f or a number of "core" galaxies 
(iMilosavhevic et aP 2002t iRavi ndranath. Ho & Filipperikol 
I2002t lGrahamll2004t iFerra rese et alj |2006l) . As first empha- 
sized bv'M ilosavlievic eTalj 62002) . computing M^ef is prob- 
lematic due to the unknown form of the galaxy's luminosity 
profile before it was modified by the binary. .GrahamI (120041) 
noted that Sersic's law provides the best global fit to the lumi- 
nosity profiles of early-type galaxies and bulges and proposed 
that mass deficits be defined in terms of the deviation of the 
inner profile from the best-fitting Sersic law . This procedure 
was followed also by Ferrares e et al.l (l2006h in their study of 
Virgo galaxies using HST/ACS data. 

Figu re 111 summari z es the results from the iGrahamI (l2004l) 
and Ferrarese e t alj (12006 ^ studies. Mass deficits from 
Graham (2004) were increased by 0.072 in the logarithm to 
correct an error in the mass-to -light ratios (A. Graham, pri- 
vate communication). iGraham (2004.) gives two estimates 
of M, for each galaxy, based on the empirical correlation 





Fig. 7. — (a) Observed mass deficits from Graham (2004) (filled circles) and Ferrarese et al. (2006) (stars). Thick, thin, dashed and dotted lines show 
Mdef/M, = 0.5, 1,2 and 4 respectively, (b) Histogram of M^ef/M,. 
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with concentrat i on par ameter (iGraham et alJl200"Tl) and on the 
iGebhardtetal! (l20()0h version of the M, — a relation. We 
recomputed SBH masses for Graham's galaxie s using the 
most curr ent version of the M , — a relation ( Ferrar ese & FordI 
l2005h . IFerrarese et alJ (|20 06) computed M, for their sample 
galaxies in the same way; we have replaced M, for three of 
their galaxies (NGC4486, NGC4374, NGC4 649) with values 
derived from detail ed kinematical modelling ( Macchetto et alJ 
Il997t lBower etal.lfr9 98: Gebhardt et al. 2001). One "core" 
galaxv from IFerrarese et al.l (12006,) - NGC 798 - was ex- 
cluded since it contains a large-scale stellar disk. 

The histogram of Mdef/M, values is shown in Figure |7j). 
There is a clear peak at Mdef/M, « 1 . 

Sersic's law - which describes the projected, or surface, 
luminosity density of early-type galaxies - implies a space 
density that varies as i{r) ^ r^"^ at small radii, where y « 
(« — 1 ) /n and n is the Sersic index. The galaxies plotted in 
Figure |2] have n > 5, hence 0.8 ^ y ^ 1.0; in other words, 
the mass deficits for these galaxies were computed under the 
assumption that the pre-existing nuclear density was a power 
law with y w 1 . This is approximately the same inner depen- 
dence as in our y = 1 Dehnen-model galaxies. While the cu- 
mulative A^-body mass deficits shown in Figure|6]and Table 2 
were based on initial models with y = 1.5, Figure |4] suggests 
that the results would have been almost identical for y = 1, at 
least during the first few stages of the merger hierarchy. 

Comparing Figures |6] and [T] we therefore conclude that 
most "core" galaxies have experienced 1 ^ 5V; ^ 3 mergers 
(i.e. 0.5 ^ Mcief /M, <y 1.5), with 9{, =2 the most common 
value. 

A few galaxies in Figure [T] have significantly larger 
mass deficits; the most extreme object is NGC 5903 with 
Mdef/M, K, 4.5. Of course the uncertaintities associated 
with the Mdef and M, values in Figure [T] may be large; 
the former due to uncertain M /L corrections, the latter due 
to various difficulties associated with SBH mass estimation 
(e.g. Maciejewski & Binney 2001; Valluri, Merritt & Em- 
sellem 2004). In the case of NGC 5903, for which M. 
was computed from the M, — a relation, we note that pub- 
lishe d velocity d ispersions for this galaxy range from 192 km 
s"' (ISmith et al . 2000) to 245 km s"^ (Davies et al. 1987); 
the corresponding range in the inferred black hole mass is 
1.36 X \(f-<,M,/MQ <,AA5 x 10** and the range inM,/^//M. 
is 1.9 ^ Mdef/M, <, 6.2. Even if a were known precisely, 
there will always be a substantial uncertainty associated with 
any M, value derived from an empirical scaling relation. 

But assuming for the moment that the numbers plot- 
ted in Figure |7] are accurate, how might the galaxies with 
Mdef/M, > 2 be explained? There are several possibilities, 
(a) These galaxies are the product of iAt 4 mergers. Two 
of the largest Mdef /M, values in Figure|7]are associated with 
M87 and M49, both extremely luminous galaxies that could 
have experienced multiple mergers. However, Figure |6] sug- 
gests that Mdef /M, values ^ 3 might be difficult to produce 
via repeated mergers, (b) The primordial density profiles in 
these galaxies were flatter than Sersic's law at small radii, 
(c) The two SBHs failed to efficiently coalesce during one or 
more of the merger events, so that a binary was present when 
a third SBH fell in. Indeed this is likely to be the case in the 
largest (lowest density) galaxies in Figure|7] which have little 
gas and extremely long timescales for loss-cone repopulation 
by two-body relaxation. Infall of a third SBH onto an uncoa- 
lesced binary is conducive to smaller values of M,, since one 



or more of the SB Hs could eventually be eject ed by the grav- 
itational slingshot (iMikkola & Valtorienll 19901) : and to larger 
values of Mdef, since multip le SBHs are more ef ficient than 
a binary at displacing stars (iMerritt et al.ll2004bl) . (The for- 
mer point is only relevant to galaxies - unlike NGC 5903 - 
for which M, has actually been measured.) The chaotic inter- 
action between three SBHs would probably also assist in the 
gravitational wave coalescence of the two most massiv e holes 
by ind ucing random changes in their relative orbit (Blae s et aP 
I2OO2I) . (d) The gravitational-wave rocket effect is believed 
capable of delivering kicks to a coalesc ing binary as large 
as ~ 200 km s~^ jp avata et al. 2004; Blanchet etal.1 l2005t 
iHerrmann et al.ll2b06t .Baker et al.,i2006) . The stellar density 
drops impulsively when the SBH is kicked out, and again 
when its orbit decays via dynamical friction. Mass deficits 
produced in this way can be as larg e as M, ( iMerritt et al.l 
2004a; Bov lan-Kolchin et al.ll2004r . (e) Stars bound to the 
infalling SBH - which were neglected in the A^-body simu- 
lations presented here - might affect Mdef, although it is not 
clear what the direction or magnitude of the change would be. 

Focussing again on the majority of galaxies in Figure |7] 
with 0.5 ^ Mdef/M, <i 1.5, we can ask whether values of 
9\[ in the range 1 < 5\t < 3 are consistent with hierarchi- 
cal models of galaxy formation. If the seeds of the cur- 
rent SBHs were present at large redshift, the ancestry of 
a bright galaxy could include dozens of mergers involv- 
ing binary SBHs dVolonteri et al.ll2003al; ISesana et al.i r2004). 
However the more relevant quantity is probably the num- 
ber of mergers since the era at which most of the gas was 
depleted, since star formation from ambient gas could re- 
generate a^ensity_£u^2_^^ a binary SBH 
(lGrahamll2004 . iHaehnelt & KauffmannI (l2002l) calculate just 
this quantity, based on semi-analytic models for galaxy merg- 
ers that include prescriptions for star- and SBH-formation. 
Their Figure 2 shows probability distributions for as a 
function of galaxy luminosity for mergers with q > 0.3. For 
galaxies Uke those in Figure |7] (My <, -21, Mb ^ -20), 
IHaehnelt & Kauffmariiil ( |2002|) find a median 9{, of 1, but 
particularly among the brightest galaxies, they find that val- 
ues of 9\[ as large as 3 or 4 are also likely. This prediction is 
consistent with Figure|2] 

Probably the biggest uncertainty in this analysis is the un- 
known behavior of the binary after it reaches a w Osndi and be- 
fore the next SBH falls in. If mergers are "dry," gas-dynamical 
torques can not be invoked to accelerate the coalescence, im- 
plying a complicated interaction between three SBHs when 
the next merger event occurs. As discussed above, the net ef- 
fect would be an increase in Mjef/M,. On the other hand, if 
gas is present in sufficient quantities to assist in coalescence, 
it may also form new stars - decreasing Mdef - and/or accrete 
onto the SBH - increasing M, ; in either case, Mdef / M, would 
be smaller than computed here. 

7. STALLED BINARY SEPARATIONS 

Equation (fT2l i gives an estimate for the stalling radius of 
a binary SBH in terms of its influence radius r^; the latter 
is defined as the radius of a sphere containing a stellar mass 
equal to twice M, = Mi +M2, after the binary has reached 
dstaii, i-C- after it has finished modifying the stellar density 
profile. One can use Equation (fT2] l to estimate binary separa- 
tions in galaxies where M, and p(r) are known, under the as- 
sumption (discussed in more detail below) that no additional 
mechanism has induced the binary to evolve beyond fl.s,a//. 
Table 3 gives the results for the seven Virgo "core" galax- 
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TABLE 3 
Virgo "Core" Galaxies 



Galaxy 
(1) 



NGC 4472 
NGC 4486 
NGC 4649 
NGC 4406 
NGC 4374 
NGC 4365 
NGC 4552 



(2) 



M. 

(3) 



(4) 



-21.8 
-21.5 
-21.3 
-21.0 
-20.8 
-20.6 
-20.3 



5.94 
35.7 
20.0 
4.54 
17.0 
4.72 
6.05 



130. (1.6) 
460. (5.7) 
230. (2.9) 
90. (1.1) 
170. (2.1) 
115. (1.4) 
73. (0.91) 



"stall 

q = 0.5 
(5) 



5.6 (0.070) 
20. (0.25) 
10. (0.13) 
4.0 (0.050) 
7.6 (0.094) 
5.0 (0.063) 
3.2 (0.040) 



"stall 

q = 0.\ 
(6) 



2.1 (0.026) 
7.6 (0.095) 
3.8 (0.047) 
1.5 (0.019) 

2.8 (0.035) 

1.9 (0.023) 

1.2 (0.015) 



astaii depend on geometry? (b) If a binary continues to harden 
below a^taii, what would be the effect on Mii^fl 
The answer to the latter question is complex. Con- 

oidor a g,alaxy containing a stalled binary, and suppose that 

Ve q ^Q®^ i^ec^^nism - gravitational encounters, time-dependent 
(7) tsyques ff0m a passing galaxy, perturbations from a third 
l^H, nt^^star formation, etc. - has the net effect of plac- 
j^^^' addJ^Sial stars on orbits that intersect the binary. These 
igtgrs wpjglje ejected, and the binary will shrink. However 
l&ere nadd5not be any net change in the mass deficit, since 
Hhl new^^l&s are first added, then subtracted, from the nuclear 
ensity .^\^eco nd -order changes in p due to the time depen- 



"352" 
733. 
776. 
590. 
832. 
533. 
757. 



Notes. - Col. (1): New General Catalog (NGC) numbers. Col. (2): Absolute 
B-band galaxy magnitudes. Col. (3): Black hole masses in 10**M,v. Col. 
(4): Black hole influence radii, defined as the radii containing a mass in stars 
equal to 2M,, in pc (arcsec). Col. (5): Binary stalling radii for q = 0.5, in pc 
(arcsec). Col. (6): Binary stalling radii for ^ = 0. 1, in pc (arcsec). Col. (7): 
Typical ejection velocity from a binary SBH with a = a^iaii (km s^'). Cols. 
(8), (9): Escape velocity (km s^'). 

ies shown in Figure [T] Influence radii were computed using 
parametric ( "core-Sersic") f i ts to t he luminosity profiles, as 
described in iFerrarese et al.l (l2006h : that paper also gives the 
algorithm which was used for converting luminosity densities 
into mass densities. When discussing real galaxies, r^^ is just 
the currently-observed influence radius r/, and henceforth the 
prime is dropped. 

Stalling radii are given in Table 3 assuming q = 0.5 and 
q ~ 0.1; fljta// for other values of q can be computed from 
Equation ( fT2] i. Table 3 also gives angular sizes of the binaries 
assuming a distance to Virgo of 16.52 Mpc. Typical sepa- 
rations between components in a stalled binary are found to 
be ~ 10 pc (0.1") (q = 0.5) and - 3 pc (0.03") {q = 0.1). 
Of course these numbers should be interpreted as upper lim- 
its; nevertheless they are large enough to suggest that binary 
SBHs might be resolvable in the brighter Virgo galaxies if 
both SBHs are luminous. 



8. DISCUSSION 

We have shown that the mass deficit produced by a binary 
SBH at the center of a spherical galaxy is "quantized" in units 
of 0.5 (Ml +M2), with only a weak dependence on binary 
mass ratio, and that M^ief /M, grows approximately linearly 
with the number of merger events. These results were com- 
pared with observed mass deficits to conclude that most bright 
elliptical galaxies have experienced 1—3 mergers since the 
era at which gas was depleted and/or star formation became 
inefficient. 

In this section, we explore some further implications of the 
A^-body results, and discuss their generality. 

8.1. Loss-Cone Refilling and Non-Spherical Geometries 

The approach adopted in this paper was motivated by the 
fact that mass deficits, or cores, are only observed at the cen- 
ters of galaxies with very long relaxation times. In such galax- 
ies, collisional repopulation of orbits depleted by the binary 
would act too slowly to significantly influence the binary's 
evolution (IValtonenlll996h . and the binary would stop evolv- 
ing once it had interacted with all stars on intersecting orbits. 
But even in the completely collisionless case, evolution of the 
binary, and its effect on the local distribution of stars, could 
be different in nonspherical (axisymmetric, triaxial, ...) ge- 
ometries. There are really two questions here: (a) How does 



dence of the gravitational potential are being ignored.) This 
argument suggests that the good correlation found here (Fig- 
ure|4]i between astaii and Mjef need not apply more generally, 
in time-dependent or non-spherical situations. 

An example of the latter is a binary embedded in a triax- 
ial galaxy which contains centrophilic (box or chaotic) or- 
bits. The mass in stars on centrophilic orbits can be extremely 
large, M,, in a triaxial galaxy, although many orbital pe- 
riods will generally be required before any given star passes 
near enough to the binary to interact with it ( Merritt & PoonI 
120041; iHollev-Bockelmann & Sigurdssonl2006). Nevertheless 
the binary need not stall (^erczik et al. 2006). Self-consistent 
calculations of the effect of a binary on the nuclear density 
profile have not been carried out in the triaxial geometry, but 
it is clear that mass deficits might be smaller than in the spher- 
ical geometry since some stars ejected by the binary will be 
on orbits with very large r/,) characteristic radii. 

The influence of different geometries on fl.s,a// and M^ef will 
be investigated in future papers. 

8.2. Mass Deficits and Core Radio Power 
An intriguing r elation wa s discovered by 



flmaverd e & Capettil (l2006h and ICapetti & Balmaverdd 
)06) between the radio and morphological properties of 
active galaxies. Radio-loud AGN - active nuclei with a large 
ratio of radio to optical or radio to X-ray luminosities - are 
uniquely associated with "core" galaxies. In fact, the slope 
of the inner brightness profile appears to be the only quantity 
that reliably predicts whether an AGN is radio-loud or 
radio-quiet. iCapetti & Balmaverdd (l2006l) also showed that 
when the radio-loud and radio-quiet galaxies are considered 
separately, there is no dependence of radio loudness on the 
nuclear density profile within either class. In other words, 
the fact that an active galaxy is morphologically a "core" 
galaxy predicts that it will be radio-loud but does not predict 
how loud. While the origin of this connection is unclear, it 
might reflect the influen ce of SBH rotation on radio power 
(IWilson & Colberti[T995h : ' 'core" galaxies have experienced 
a recent merger, and coalescence of the two SBHs resulted in 
a rapidly-spinning remnant. 

Two results from the A^-body work presented here may 
be relevant to the Capetti & Balmaverde correlation. First, 
as shown above, the mass deficit is a weak function of the 
mass ratio of the binary that produced it, i.e., the shape 
of the nuclear density profile is poorly correlated with the 
mass of the infalling hole. Since the degree of spin-up 
is also weakly correlated \y i th M7/M1 (IW ilson & Colberj 
[T99I iMerritt & EkersI I200I iHughes & Blandford. ,2003;) - 
even very small infalUng holes can spin up the larger hole 
to near-maximal spins - it follows that SBH rotation should 
be poorly correlated with mass deficit. This may explain the 
weak dependence of radio loudness on profile slope found by 
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ICapetti & Balmaverdd (l2006h for the core galaxies. 

Second, as argued above, core galaxies may contain uncoa- 
lesced binaries. Perhaps the presence of a second SBH is the 
main factor that determines radio loudness. 

8.3. High-Velocity Stars 

Gravitational slingshot ejections by a binary SBH produce 
a population of high- velocity stars with trajectories directed 
away from the center of the galaxy. A few hyper- velocity stars 
have been detected in the halo of the Milky Wa y jBrown et al.l 
l2005t iHirsch et all l2005t iBrown et all l2006l) . and a binary 
SBH at the Galactic center is a possible model for their ori- 
gin. The short nuclear relaxation time would accelerate the 
hardening of a massive binary by e nsuring that stars wer e 
scattered into its sphere of influence dYu & Tremainell2003h . 
In this model, stars could have been ejected during the late 
stages of the binary's evolution, when a ^ a/,, at high enough 
velocities and large enough numbers to explain the observed 
hyper- velocity stars. 

Binary SBHs have also been invoked to explain other popu- 
lati ons, e.g. intergalactic p lanetary nebulae in the Virgo clus- 
ter dArnaboldi et al.ll2004l) . If every Virgo galaxy harbored a 
binary SBH with mass ratio q^O.l which evolved, via sling- 
shot ejection of stars, to the gravitational -radiation regime, the 
total mass ejected would have been ~ 2% of the l uminous 
mass of the cluster ( Hollev-Bockelmann et alj[2005h . But at 
least in the brighter Virgo galaxies - including the galaxies 
in Table 3 - two-body relaxation times are much too long for 
collisional resupply of binary loss cones, and there is no com- 
pelling reason to assume that the binary SBHs in these galax- 
ies would have continued interacting with stars beyond a^ta//. 
At these separations, velocities between the two components 
of the binary are relatively modest, implying a much lower 
probability of ejection of stars at velocities large enough to 
escape the galaxy. 

The mean specific energy chang e of a star that intera cts with 
a hard binary is ~ 3G^/2fl (Hifls 1 983HOuinlan|[T996l) . Com- 
bining this with Equation ( fT2l i for Osraii gives the typical ejec- 
tion velocity from a stalled binary: 



i4.0 



1/2 



1/2 



rh 



value of M,. Ejection velocities are seen to be much smaller 
than \esc in all cases, implying that essentially no stars would 
be ejected into the intracluster medium. 

8.4. Dark-Matter Cores 

The specific energy change of a dark matter particle inter- 
acting with a binary SBH is identical to that of a star. At least 
in the case of bright elliptical galaxies like those in Table 3, it 
is unlikely that dark matter was ever a dominant component 
near the center or that it significantly influenced the evolu- 
tion of the binary. But ejection of dark matter particles by a 
massive binary would produce a core in the dark matter distri- 
bution similar in size to the luminous-matter core, rcore ~ Oi- 
A more definite statement about the variation with radius of 
Pdm/P* near the center of a galaxy containing a binary SBH 
is probably impossible to make without A^-body simulations 
that contain all three components. However Table 3 suggests 
dark matter core radii of hundreds of parsecs in bright ellipti- 
cal galaxies. 

It follows that the rates of self-interaction of su persymmet- 
ric pa rticles at the centers of galaxies like M87 (Baltz et alj 
I2OOOI) would be much lower than computed under the assump- 
tion th at the dark matter stiU retains its "prim ordial" density 
profile (iNavarro et al.lll996HMoore et al.lll998 1. This fact has 
implications for so-called "indirect" dark matter searches, in 
which inferences are drawn about the properties of particle 
dark matt er based on measurem ents of its self-annihilation by- 
products (iBertone et al.l 120041) . One of the proposed search 
strategies is to identify a component of the diffuse gamma- 
ray background that is gen erated by dark matter annhila 
tions in halos at all redshifts (' Ullio et al.ir2002 ': Tavlor & Sill 
2003). Calcula tions of the background flux (Ando 2005L 
lElsasser & Mann heim 2005) have so far always assumed that 
the dark matter distribution does not change with time. Since 
the annihilation flux in a smooth dark matter halo is dom- 
inated by the the center, these calculations may substantially 
over-estimate the contribution of galactic halos to the gamma- 
ray background. 



10 pc, 
(15) 

independent of mass ratio. Of course this is an upper limit in 
the sense that ejection velocities are lower when f < t stall- Ta- 
ble 3 gives Vgj for the Vu'go "core" galaxies and also v^y^, the 
escape velocity, defined as -^Z— 2<I>(fl^„//) with 4>(r) the grav- 
itational potential including the contribution from the binary, 
modeUed as a point with mass equal to the currently-observed 
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